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Abstract 



The anisotropy of the azimuthal distributions of charged particles produced in 
^SfvfN = 2.76 TeV PbPb collisions is studied with the CMS experiment at the LHC. 
The elliptic anisotropy parameter, vi, defined as the second coefficient in a Fourier 
expansion of the particle invariant yields, is extracted using the event-plane method, 
two- and four-particle cumulants, and Lee-Yang zeros. The anisotropy is presented as 
a function of transverse momentum (px)/ pseudorapidity (//) over a broad kinematic 
range: 0.3 < px < 20GeV/c, \\]\ < 2.4, and in 12 classes of collision centrality from 
to 80%. The results are compared to those obtained at lower center-of-mass energies, 
and various scaling behaviors are examined. When scaled by the geometric eccentric- 
ity of the collision zone, the elliptic anisotropy is found to obey a universal scaling 
with the transverse particle density for different collision systems and center-of-mass 
energies. 
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1 Introduction 

The azimuthal anisotropy of emitted charged particles is an important feature of the hot, dense 
medium produced in heavy-ion collisions, and has contributed to the suggestion of a strongly 
coupled quark-gluon plasma (sQGP) being produced in nucleus-nucleus collisions at RHIC [I]- 
|5l. In noncentral collisions, the beam direction and the impact parameter vector define a reac- 
tion plane for each event. If the nucleon density within the nuclei is continuous, the initial 
nuclear overlap region is spatially asymmetric with an "almond-like" shape. In this approxi- 
mation, the impact parameter determines uniquely the initial geometry of the collision, as il- 
lustrated in Fig.[lj In a more realistic description, where the position of the individual nucleons 
that participate in inelastic interactions is considered, the overlap region has a more irregular 
shape and the event-by-event orientation of the almond fluctuates around the reaction plane. 
Experimentally, the azimuthal distribution of the particles detected in the final state can be used 
to determine the "event plane" that contains both the beam direction and the azimuthal direc- 
tion of maximum particle density. Strong rescattering of the partons in the initial state may lead 
to local thermal equilibrium and the build up of anisotropic pressure gradients, which drive a 
collective anisotropic expansion. The acceleration is greatest in the direction of the largest pres- 
sure gradient, i.e., along the short axis of the almond. This results in an anisotropic azimuthal 
distribution of the final-state hadrons. The anisotropy is quantified in terms of a Fourier ex- 
pansion of the observed particle yields relative to the event-by-event orientation of the event 
plane IH: 



E— — = — — = 1+ y 2vn(pT,y) cos \n (cp -1")] , (1) 

d3p pxdpxdi/d^ iK^id-pidy \ 'nt^i'VJ i \y )\ y \! 





Figure 1: A schematic diagram of a noncentral nucleus-nucleus collision viewed in the plane or- 
thogonal to the beam. The azimuthal angle (p, the impact parameter vector b, and the reaction- 
plane angle Yr are shown. The event-plane angle Y, with respect to which the flow is mea- 
sured, fluctuates around the reaction-plane angle. 

where q), E, y, and px are the particle's azimuthal angle, energy, rapidity, and transverse mo- 
mentum, respectively, and Y is the event-plane angle. The second coefficient of the expansion. 
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often referred to as the "elliptic flow" strength, carries information about the early collision dy- 
namics IZllHI- The coefficients in the Fourier expansion may depend on pj, rapidity, and impact 
parameter. Typically, the measurements are obtained for a particular class of events based on 
the centrality of the collisions, defined as a fraction of the total inelastic nucleus-nucleus cross 
section, with 0% denoting the most central collisions with impact parameter b = 0, and 100% 
the most peripheral collisions. Expressions similar to Eq. ([T]) can be written for the Fourier 
expansion of the yield integrated over pj or rapidity: 

|^ = i;f (l + L^^-McosM.-T)]), (2, 

and 



/\ =^ — =^ l+r2z;„(pT)cos[n(^>-Y)l . (3) 

The Fourier coefficients Vn{y) and Vnipj) can be obtained by directly analyzing the yields in- 
tegrated over Pj or rapidity, or from the coefficients of the triple-differential invariant yield 
in Eq. ([T||, v„{pj,y), by performing the following yield-weighted average over the ranges of 
transverse momentum Apj and rapidity Ay, from which the particles are taken: 



vn{y) 



Pt dpT dy 



(4) 



or 



J Ay ■1' Pt dpTdy 

The coefficients of the Fourier expansion of the particles' invariant yield integrated over a broad 
rapidity and pj window, are often referred to as "integrated flow": 



dN 
d^ 

where 




X]2t;„cos[n((p-Y)] j , (6) 



_ Ly /apt Pt dpT^j^^VnjpT, y) 

JAydy J Apt PTdpT^4^ 

In obtaining the Fourier coefficients, the absolute normalization in the particle yields is not 
important, as long as the particle detection efficiency is constant over the chosen transverse 
momentum and rapidity range. However, if the efficiency varies, the appropriate efficiency 
corrections need to be applied. This often leads to a two-step procedure in which first v„{pj,y) 
is obtained in a narrow phase-space window where the efficiency is constant, and then a yield- 
weighted average is performed using Eqs. (Q, (|5|, or folding in the efficiency-corrected 
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particle spectra. When the particle mass is not determined in the measurement, the pseudora- 
pidity f] = — ln[tan(0/2)], with 9 being the polar angle, is used instead of the rapidity. 

Elliptic flow has been measured at the AGS, SPS, and RHIC. A notable feature in these mea- 
surements is that the elliptic flow measured as a function of transverse momentum, V2{pt), 
increases with the nucleon-nucleon center-of-mass energy (a/snn) up to about 22GeV [,9l[T0 | . 
and then saturates at a value compatible with predictions from ideal hydrodynamics lim[T2| . 
The most extensive experimental studies have been performed at the highest RHIC energy 
of \/snn = 200 GeV in AuAu collisions |1h3J. First results from PbPb collisions at i/s^n = 
2.76 TeV from the Large Hadron Collider (LHC) IITSllMll indicate that there is little or no change 
in the transverse momentum dependence of the elliptic flow measured at the highest RHIC en- 
ergy and the LHC, despite the approximately 14-fold increase in the center-of-mass energy. Re- 
cent theoretical studies of elliptic flow have focused on quantifying the ratio of the shear viscos- 
ity to the entropy density of the produced medium assuming viscous hydrodynamics |T5l [T6| , 
and taking into account a variety of possible initial conditions. 

Based on experimental results and the corresponding theoretical descriptions of the data, the 
underlying physics processes that generate the elliptic anisotropy are thought to vary for dif- 
ferent kinematic regions: 

• Elliptic flow in the bulk system 

Hadrons produced in soft processes, carrying low transverse momentum (pj < 
2 GeV/c for mesons, and pi ^ 3 GeV/c for baryons [T7L1201), exhibit azimuthal anisotropies 
that can be attributed to collective flow driven by the pressure gradients in the sys- 
tem. The description of elliptic flow is amenable to hydrod}mamic calculations |(T2l 
[I5II2TH22I. Comparisons to theory indicate that the flow is primarily generated dur- 
ing an early stage of the system evolution. 

• Recombination region 

At intermediate transverse particle momentum (2 < px ^ 4GeV/c), the RHIC data 
show that the elliptic anisotropies for various identified hadron species approxi- 
mately follow a common behavior when both the ^^2(^1) value and the pj of the 
particle are divided by the number of valence quarks in the hadron lllTl [181 123^. This 
behavior is successfully reproduced by models invoking quark recombination as the 
dominant hadronization mechanism in this momentum range Il24l - |26l| . The quark- 
number scaling of the elliptic flow has been interpreted as evidence that quark de- 
grees of freedom dominate in the early stages of heavy-ion collisions, when the col- 
lective flow develops | |27| . The quark recombination may involve both thermally 
produced quarks and quarks originating in jet fragmentation. Therefore, the elliptic 
anisotropy in the recombination region results from an interplay between the bulk 
flow of the system and the azimuthal anisotropy in hadron production induced by 
jet quenching. 

• Jet-fragmentation region 

At intermediate and high transverse momentum (pj > 3GeV/c), where fragments 
from increasingly harder partons begin to contribute to the particle spectra, anisotropy 
in the azimuthal distributions may be generated from the stronger jet quenching in 
the direction of the long axis of the almond-shaped reaction zone L28r.31]. It is ex- 
pected that this mechanism will dominate the elliptic anisotropy of hadrons with 
momenta in excess of 8 GeV/c. Thus, measurements extending beyond this pj range 
carry information about the path-length dependence of energy loss in the produced 
medium. 
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2 Experimental method 



The measurements presented in this paper span the transverse momentum range of 0.3 < 
Pt < 20 GeV/ c and provide the basis for comparisons to theoretical descriptions of the bulk 
properties of the system and the quenching of jets. Such comparisons may give insight in 
determining the transport properties of sQGP, namely: the shear viscosity, and the opacity of 
the plasma. The theoretical interpretation of the elliptic anisotropy in the recombination region 
requires identified-hadron measurements that are not included in this analysis. 

In ideal hydrodynamics, the integrated elliptic flow is directly proportional to the initial spa- 
tial eccentricity of the overlap zone [32|. There are many factors that can change this behavior, 
including viscosity in the sQGP and the hadronic stages of the system evolution, incomplete 
thermalization in peripheral collisions, and variations in the equation-of-state. Event-by-event 
fluctuations in the overlap zone l!33ti38| could also influence the experimental results, depend- 
ing on the method that is used to extract the V2 signal. Invoking multiple methods with dif- 
ferent sensitivities to the initial-state fluctuations is important for disentangling the variety of 
physics processes that affect the centrality dependence of the elliptic flow. Comparisons to re- 
sults from lower-energy measurements and explorations of empirical scaling behaviors could 
provide additional insights into the nature of the matter produced in high-energy heavy-ion 
collisions. 

The pseudorapidity dependence of the elliptic flow, V2{f]), provides information on the initial 
state and the early-time development of the collision, constraining the theoretical models and 
the description of the longitudinal dynamics in the collisions 11391 |40I . Longitudinal scaling in 
V2{f]) extending over several units of pseudorapidity (extended longitudinal scaling) has been 
reported at RHIC [41J for a broad range of collision energies (^Jsmn = 19.6-200 GeV). Studies 
of the evolution of ^2(1]) from RHIC to LHC energies may have implications for the unified 
description of sQGP in both energy regimes. 

The paper is organized as follows: Section |2] gives the experimental details of the measure- 
ment including the Compact Muon Solenoid (CMS) detector, triggering and event selection, 
centrality determination, Glauber-model calculations, reconstruction of the charged-particle 
transverse momentum spectra, methods of measuring the elliptic anisotropy, and the studies 
of the systematic uncertainties. In Section |3] we present results of V2 as a function of trans- 
verse momentum, centrality, and pseudorapidity, and the measurement of the charged-particle 
transverse momentum spectra. The elliptic flow results are obtained with the event-plane 
method |6|, two- and four-particle cumulants [42| , and Lee-Yang zeros Il43l |44| . The analy- 
sis is performed in 12 classes of collision centrality covering the most central 80% of inelastic 
collisions. We study the eccentricity scaling of V2, and investigate the differences in the results 
obtained from different methods, taking into account their sensitivity to initial-state fluctua- 
tions. Section |4] is devoted to detailed comparisons with results obtained by other experiments 
at lower energies and the exploration of different scaling behaviors of the elliptic flow. The 
results of our studies are summarized in Section [S] 

2 Experimental method 
2.1 CMS detector 

The central feature of the CMS apparatus is a superconducting solenoid of 6 m internal di- 
ameter, providing a 3.8 T field. Within the field volume are a silicon tracker, a crystal elec- 
tromagnetic calorimeter, and a brass /scintillator hadron calorimeter. Muons are measured in 
gas-ionization detectors embedded in the steel return yoke. In addition to these detectors, 
CMS has extensive forward calorimetry. The inner tracker measures charged particles within 
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the pseudorapidity range < 2.5, and consists of silicon pixel and silicon strip detector mod- 
ules. The beam scintillation counters (BSC) are a series of scintillator tiles which are sensitive 
to almost the full PbPb interaction cross section. These tiles are placed along the beam line 
at a distance of ±10.9 m and ±14.4 m from the interaction point, and can be used to provide 
minimum-bias triggers. The forward hadron calorimeter (HF) consists of a steel absorber struc- 
ture that is composed of grooved plates with quartz fibers inserted into these grooves. The HF 
calorimeters have a cylindrical shape and are placed at a distance of 11.2 m from the interaction 
point, covering the pseudorapidity range oi 2.9 < \ri\ < 5.2. A more detailed description of the 
CMS detector can be found elsewhere ||45l . 

The detector coordinate system has the origin centered at the nominal collision point inside 
the experiment, with the y axis pointing vertically upward, the x axis pointing radially inward 
towards the center of the LHC ring, and the z axis pointing along the counterclockwise beam 
direction. 

2.2 Event selection 

The measurements presented are performed by analyzing PbPb collision events recorded by 
the CMS detector in 2010. From these data, the minimum-bias event sample is collected using 
coincidences between the trigger signals from both the ±z and — z sides of either the BSC or 
HF. The minimum-bias trigger used for this analysis is required to be in coincidence with the 
presence of both colliding ion bunches in the interaction region. This requirement suppresses 
noncollision-related noise, cosmic rays, and beam backgrounds (beam halo and beam-gas colli- 
sions). The total hadronic collision rate varied between 1 and 210 Hz, depending on the number 
of colliding bunches and the bunch intensity. 

In order to obtain a pure sample of inelastic hadronic collisions, several offline selections are ap- 
plied to the triggered event sample. These selections remove contamination from noncollision 
beam backgrounds and from ultraperipheral collisions (UPC) that lead to an electromagnetic 
breakup of one or both of the Pb nuclei [46 J. First, beam-halo events are vetoed based on the 
BSC timing. Then, to remove UPC and beam-gas events, an offline HF coincidence of at least 
three towers on each side of the interaction point is required, with a total deposited energy of 
at least 3 GeV in each tower. A reconstructed primary vertex made of at least two tracks and 
consistent with the nominal interaction point position is required. To further reject beam-gas 
and beam-halo events, the pixel clusters are required to have a length along the beam direc- 
tion compatible with being produced by particles originating from the primary vertex, as for 
the study in [47]. Additionally, a small number of noisy events with uncharacteristic hadron 
calorimeter responses are removed. 

The acceptance of the silicon tracker for \rj\ > 0.8 is found to vary significantly with respect to 
the longitudinal position of the collision point relative to the geometric center of the detector. 
This event-by-event variance in the tracking efficiency contributes to a systematic bias of the 
elliptic flow measurements at forward pseudorapidity. In order to remove this bias, events in 
this analysis are required to have a longitudinal vertex position within 10 cm of the geometric 
center of the detector. After all selections, 22.6 million minimum-bias events, corresponding to 
an integrated luminosity of approximately 3 }ih^^, remain in the final sample. 

2.3 Centrality determination and Glauber model calculations 

In this analysis, the observable used to determine centrality is the total energy in both HF 
calorimeters. The distribution of this total energy is used to divide the event sample into 40 cen- 
trality bins, each representing 2.5% of the total nucleus-nucleus interaction cross section. The 
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events are then regrouped to form 12 centrality classes used in the analysis: 0-5% (most cen- 
tral), 5-10%, 10-15%, 15-20%, 20-25%, 25-30%, 30-35%, 35^0%, 40-50%, 50-60%, 60-70%, and 
70-80% (see Table [l]|. Using Monte Carlo (MC) simulations, it is estimated that the minimum- 
bias trigger and the event selections include (97 ± 3)% of the total inelastic cross section. For 
the events included in this analysis (0-80% centrality), the trigger is fully efficient. 

For each group of events that comprises a centrality class, we evaluate a set of quantities that 
characterize the initial geometry of the collisions using a MC Glauber model. The Glauber 
model is a multiple-collision model that treats a nucleus-nucleus collision as an independent 
sequence of nucleon-nucleon collisions (see Ref. [48 1 and references therein). The nucleons 
that participate in inelastic interactions are called "participants". A schematic view of a PbPb 
collision with an impact parameter b = 6 fm, as obtained from the Glauber model, is shown in 
Fig. |2j The direction and the magnitude of the impact parameter vector and the corresponding 
reaction-plane angle Yr, are the same as in Fig. [l] However, the initial interaction zone as 
determined by the spatial distribution of the participants (filled circles) is no longer regular in 
shape and is not necessarily symmetric with respect to the reaction plane. In each event one 
can evaluate the variances cr^ and a^, and the covariance axy = {xy) — [x) [y) of the participant 
distributions projected on the x and y axes. One can then find a frame x'-y' that minimizes 
a^i, and define a "participant plane" using the beam direction and the x' axis 1341 149| . In this 
frame, the covariance (J^^'y' of the participant spatial distribution vanishes. To characterize the 
geometry of the initial state of the collision, we define Il34ll49l the eccentricity of the participant 
zone epart/ its cumulant moments e{2} and e{4}, and the transverse overlap area of the two 
nuclei S: 



-part 



e{2}2 ^ {el,,,), 



e{4}4 ^ 2(e2,,,)2 - (4,,,), and 



,4 
part/ 



(8) 

(9) 
(10) 



na^, ay, = nJala^-aly. (11) 



x^y ^xy 

In Eqs. (|9]) and ( 10 1, the average is taken over many events from the same centrality interval. 

To implement the Glauber model for PbPb collisions at y/swti = 2.76 TeV, we utilize the founda- 
tion of a published Glauber MC software package TGlauberMC [50J, which was developed 
for the PHOBOS Collaboration at RHIC. Standard parameters of the Woods-Saxon function 
used for modeling the distribution of nucleons in the Pb nuclei are taken from Ref. [51J. The 
nucleon-nucleon inelastic cross section, which is used to determine how close the nucleon tra- 
jectories need to be in order for an interaction to occur, is taken as 64 ± 5 mb, based on a fit 
of the existing data for total and elastic cross sections in proton-proton and proton-antiproton 
collisions II52I . The uncertainties in the parameters involved in these calculations contribute to 
the systematic uncertainty in Npart/ S, and Cpart for a given centrality bin. The connection be- 
tween the experimentally defined centrality classes using the HF energy distribution and Npart 
from the Glauber model is obtained |[53i| from fully simulated and reconstructed MC events 
generated with the AMPT ||54| event generator. The calculated Glauber model variables for 
each centrality class are shown in Table [l] 
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Figure 2: (Color online) A schematic view of a PbPb collision with an impact parameter b = 
6 fm as obtained from the Glauber model. The nucleons that participate in inelastic interactions 
are marked with filled circles. The x and y coordinates represent the laboratory frame, while x' 
and y' represent the frame that is aligned with the axes of the ellipse in the participant zone. The 
participant eccentricity gpart and the standard deviations of the participant spatial distribution 
CTyi and CTx' from which the transverse overlap area of the two nuclei is calculated are also shown. 
The angle Yr denotes the orientation of the reaction plane. 



2.4 Reconstruction of the charged-particle transverse momentum distributions 
and the mean transverse momentum 

To determine the transverse momentum distributions of the charged particles produced in the 
collisions, we first need to reconstruct the particles' trajectories ("tracks") through the 3.8 T 
solenoidal magnetic field. The tracks are reconstructed by starting with a "seed" compris- 
ing two or three reconstructed signals ("hits") in the inner layers of the silicon strip and pixel 
detectors that are compatible with a helical trajectory of some minimum pj and a selected re- 
gion around the reconstructed primary vertex or nominal interaction point. This seed is then 
propagated outward through subsequent layers using a combinatorial Kalman-filter algorithm. 
Tracking is generally performed in multiple iterations, varying the layers used in the seeding 
and the parameters used in the pattern recognition, and removing duplicate tracks between 
iterations. This algorithm is described in detail in Ref. [55J. The algorithm used in most of the 
CMS proton-proton analyses, as well as the tracking detector performance for the 2010 run, are 
described in Ref. [56]. 

The six-iteration process used in proton-proton collisions is computationally not feasible in the 
high-multiplicity environment of very central PbPb collisions. In place of this, a simple two- 
iteration process is used. The first iteration builds seeds from hits in some combination of three 
layers in the barrel and endcap pixel detectors compatible with a trajectory of px > 0.9GeV/c 
and a distance of closest approach to the reconstructed vertex of no more than 0.1 cm in the 
transverse plane and 0.2 cm longitudinally. These tracks are then filtered using selection crite- 
ria based on a minimum number of reconstructed hits, vertex compatibility along the longitu- 
dinal direction and in the transverse plane, and low relative uncertainty on the reconstructed 
momentum. 
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Table 1: For the centrality bins used in the analysis, the average values of the number of partic- 
ipating nucleons, transverse overlap area of the two nuclei, participant eccentricity, and cumu- 
lant moments of the participant eccentricity, along with their systematic uncertainties from the 
Glauber model. 
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0-5% 


381±2 


29.4±1.2 


0, 


.074±0.003 


0, 


.084±0. 


,003 


0. 


,053±0, 


.002 


5-10% 


329±3 


26.6±1.1 


0, 


.111±0.005 


0, 


.128±0. 


,005 


0. 


,060±0, 


.003 


10-15% 


283±3 


24.0±1.0 


0, 


.154±0.007 


0, 


.175±0. 


,007 


0. 


,122±0, 


.005 


15-20% 


240±3 


21.6±1.0 


0, 


.198±0.009 


0, 


.219±0. 


,009 


0. 


,171 ±0, 


.007 


20-25% 


204±3 


19.5±0.9 


0, 


.238±0.009 


0, 


.262±0. 


,010 


0. 


,214±0, 


.008 


25-30% 


171±3 


17.5±0.8 


0, 


.276±0.010 


0, 


.302±0. 


,012 


0. 


,253±0, 


.010 


30-35% 


143±3 


15.7±0.8 


0, 


.312±0.011 


0, 


.339±0. 


,012 


0. 


,288±0, 


.010 


35-40% 


118±3 


14.1±0.7 


0, 


.346±0.010 


0, 


.375±0. 


Oil 


0. 


,322±0, 


.009 


40-50% 


86.2±2.8 


12.0±0.6 


0, 


.395±0.010 


0, 


.429±0. 


Oil 


0. 


,370 ±0, 


.010 


50-60% 


53.5±2.5 


9.4±0.5 


0, 


.465±0.008 


0, 


.501±0. 


,009 


0. 


,437±0, 


.007 


60-70% 


30.5±1.8 


7.1±0.4 


0, 


.543±0.011 


0, 


.581±0. 


,012 


0. 


,514±0, 


.010 


70-80% 


15.7±1.1 


4.8±0.3 


0, 


.630±0.016 


0, 


.662±0. 


,017 


0. 


,598±0, 


.015 



In the second iteration, seeding is also performed using three layers of the pixel detector, but the 
minimum transverse momentum requirement is relaxed to pj > 0.2 GeV/c. These tracks are not 
propagated through the silicon-strip detector, but simply refitted using the transverse position 
of the beam spot as an additional constraint. These pixel-only tracks are then filtered using 
selection criteria of vertex compatibility along the longitudinal axis and statistical goodness of 
fit. 

The tracks from both collections are checked for duplicate tracks using the number of hits in 
common between the two tracks, and duplicates are removed giving preference to the first- 
iteration tracks. The tracking algorithm may sometimes misidentify tracks by combining sil- 
icon detector signals that do not originate from the same charged particle. It is important to 
keep the proportion of misidentified tracks (referred to as "fake tracks"), or fake rate, as low as 
possible. To create the final collection, first-iteration tracks with pj > 1.5GeV/c are combined 
with second-iteration pixel-only tracks with pj < 1.8GeV/c. These px limits were chosen to 
exclude kinematic regions where a given iteration has a high fake rate. The efficiency and fake 
rate of this modified tracking collection are found both by using a full MC simulation of PbPb 
collisions based on the hydjet event generator II57I and by embedding simulated charged pi- 
ons into PbPb data events. The efficiency, fake rate, and momentum resolution of this final 
collection determined by the simulated events from the hydjet event generator are shown in 
Fig. |3] for events of five different centrality classes. The abrupt change in efficiency, fake rate, 
and momentum resolution seen in the figure occurs at the transverse momentum where the 
two track collections are merged. 

The charged-particle transverse momentum distributions are corrected for the loss of accep- 
tance, efficiency, and the contributions from fake tracks. Each detected track is weighted by a 
factor Wtr according to its centrality, transverse momentum, and pseudorapidity: 

1 - f 

iffr (centrality, px/ '/) = — (12) 

Here / is the fraction of fake tracks that do not correspond to any charged particle, and e is the 
absolute tracking efficiency accounting for both geometric detector acceptance and algorith- 
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Figure 3: (Color online) Efficiency (top), fake rate (middle), and momentum resolution (bottom) 
of charged tracks obtained from hydjet simulated events in four pseudorapidity regions: \ti\ < 
0.8, 0.8 < \r]\ < 1.6, 1.6 < \t]\ < 2.0, and 2.0 < < 2.4 displayed from left to right, and for the 
five centrality classes given in the legend. 
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mic tracking efficiency. The proportion of multiply reconstructed particles and reconstructed 
particles from secondary decays is negligible and is not included in the correction factor. The 
fully corrected transverse momentum distributions are measured in 12 centrality classes over 



the pseudorapidity range \rj\ < 2.4 in bins of Arj = 0.4, as discussed in Section 3.4 These dis- 
tributions are used in obtaining integrated V2 values. We also study the evolution of (px) with 
pseudorapidity and centrality (Section [3!4] |, and center-of-mass energy (Section |4]|. 

To evaluate {pj), the spectra need to be extrapolated down to pj = OGeV/c. The extrapolation 
is performed using a Tsallis distribution 



dp'^ 27rpr p arjapr \ nT J 

where Ej = + pj — m, and m is taken to be the charged pion mass. The measured spectra 
are fitted in the range 0.3 < pj < 3.0 GeV/c , and the fit parameters C, n, and T are determined. 
The mean transverse momentum is then evaluated using the fit function in the extrapolation 
region of < px < 0.3GeV/c, and the data from the range 0.3 < px < 6.0GeV/c. This method 
has been previously applied in CMS [47J in the measurement of (px) in pp collisions at ^/s = 
7TeV. 

2.5 Methods for measuring the anisotropy parameter V2 

Anisotropic flow aims to measure the azimuthal correlations of the particles produced in heavy- 
ion collisions in relation to the initial geometry of the collisions. Originally, the flow was de- 
fined as a correlation of the particle emission angles with the reaction plane HHZl. More recently, 
it was recognized I|34ll49l that the initial geometry is better characterized by the positions of the 
individual nucleons that participate in inelastic interactions, and thus define a participant plane 
that fluctuates around the reaction plane on an event-by-event basis. Neither the reaction plane 
nor the participant plane are directly measurable experimentally. Instead, there are several ex- 
perimental methods that have been developed to evaluate the anisotropic flow based on the 
final-state particle distributions. In the present analysis, we use the event-plane method, two- 
particle and four-particle cumulants p2l, and the Lee-Yang zeros method Il43l |44| the last of 
which is based on correlations involving all the particles in the event. 

The anisotropic flow measurements are affected by fluctuations that come from several sources. 
Statistical fluctuations arise due to the fact that a finite number of particles is used to determine 
a reference plane for the flow and the multiplicity fluctuations within the chosen centrality 
interval. The effect of these fluctuations is to reduce the measured flow signal and is largely 
compensated for by the resolution corrections described below. Any remaining effects from 
statistical fluctuations on our measurements are included in the systematic uncertainties (see 



Section 2.7 1. Another more important source of fluctuations comes from the event-by-event 
fluctuations in the participant eccentricity that are present even at fixed Npart- These dynamical 
fluctuations have been shown ||34| -|38| to affect the various methods for estimating the flow 
differently, since each of them is based on a different moment of the final-particle momentum 
distribution. For example, the two-particle cumulant method measures an r.m.s. value of the 
flow that is higher than the mean value. Conversely, the four-particle cumulant and other 
multiparticle correlation methods return a value that is lower than the mean value. For the 
event-plane method, the results vary between the mean and the r.m.s. value depending on the 
event-plane resolution, which varies in each centrality interval. 

In addition, there exist other sources of correlations in azimuth, such as those from resonance 
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decays, jets, and Bose-Einstein correlations between identical particles. These correlations, 
which are not related to the participant plane, are called nonflow correlations. The various 
methods proposed to estimate the magnitude of anisotropic flow have different sensitivities to 
the nonflow correlations, thus allowing systematic checks on the flow measurements. The mul- 
tiparticle correlation methods are least affected by nonflow correlations, but they do not work 
reliably when either the flow anisotropy (^2) or the multiplicity in the selected phase-space 
window is small. This happens in the most central and in the most peripheral events flZ^.'lSl. 
Thus, the two-particle cumulant and the event-plane methods provide an extended centrality 
range, albeit with a larger nonflow contribution. 

All four methods used here have been extensively studied and applied in different experi- 
ments. Thus, we limit our description in the following subsections only to the features that are 
specific for our implementations of these methods. 

2.5.1 Event-plane method 

The event-plane method estimates the magnitude of anisotropic flow by reconstructing an 
"event-plane" containing both the beam direction and the direction of the maximal flow de- 
termined from the azimuthal distributions of the final-state particles. Under the assumption 
that the flow is driven by the initial-state asymmetry in the nuclear overlap zone, and that 
there are no other sources of azimuthal correlations in the final-state particles, the event-plane 
is expected to coincide with the participant plane [34] defined in Fig.|2] Recent theoretical calcu- 
lations [38l l6TI - l63l confirm that the event plane and the participant plane are strongly correlated 
event by event. Since the event plane is determined using a finite number of particles, and de- 
tected with finite angular resolution, the measured event-plane angle fluctuates about its true 
value. As a result, the observed particle azimuthal anisotropy with respect to the event-plane 
is smeared compared to its true value. The true elliptic flow coefficient V2 in the event-plane 
method is evaluated by dividing the observed z?*® value by a resolution correction factor, R, 
which accounts for the event-plane resolution. 

To determine the event-plane resolution correction, a technique that sorts the particles from 
each event into three subevents based on their pseudorapidity values ||6| is used. For subevents 
A, B, C in three different pseudorapidity windows, the event-plane resolution correction factor 
Ra for subevent A is found as: 



where W^, Y^, and Y*- are the event-plane angles determined from the corresponding subevents, 
and the average is over all events in a selected centrality class used in the V2 analysis. 

In our implementation of the method, the event-plane angle determined for the subevent fur- 
thest in Tj from the track being used in the elliptic flow analysis is used, and the correspond- 
ing resolution correction is employed. This selection minimizes the contributions of auto- 
correlations and other nonflow effects that arise if the particles used in the event-plane de- 
termination and those used in the flow analysis are close in phase space. 

To achieve the largest pseudorapidity gap possible, two event planes are defined, with calori- 
meter data covering the pseudorapidity ranges of — 5 < ?/ < — 3 and 3 < < 5, labeled "HF-" 
and "HF+", respectively. These pseudorapidity ranges are primarily within the coverage of the 
HF calorimeters. A third event plane, found using charged particles detected in the tracker 
in the pseudorapidity range —0.8 < t] < 0.8, is also defined and used in the three-subevent 




(14) 
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technique for determining the resolution corrections for HF- and HF+. The resulting resolu- 
tion corrections are presented in Fig. |4j Particles detected in the tracker with ?] > are then 
correlated with the HF- event plane, and those with tj < with the HF+ plane. In this man- 
ner, the minimum pseudorapidity gap between particles used in the event-plane determination 
and those for which the V2 signal is measured is 3 units. A two-subevent technique based on 
HF- and HF+, with a resolution parameter defined as Ra/b = V (cos[2(Y"^ — Y^)]), is also im- 
plemented and used for systematic studies. The values of Ra/b are shown for comparison in 
Fig.g 

A standard flattening procedure [64J using a Fourier decomposition of the distribution of the 
event-plane angles to order is used to shift the event-by-event plane angle to correct for 
asymmetries in the event-plane distribution that arise from the detector acceptance and other 
instrumental effects. Although most of these effects are already accounted for with a correction 
involving just the first four coefficients in the expansion, the larger order was used for data 
quality monitoring purposes. For each centrality class in the analysis, the flattening parame- 
ters are calculated by grouping events according to the longitudinal location of their primary 
collision vertex in 5 cm-wide bins. 
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Figure 4: (Color online) Event-plane resolution correction factors as a function of centrality for 
the two event-planes (HF- and HF+) used in determining the elliptic anisotropy parameter V2- 
The corrections determined with the three-subevent method used in the analysis are shown as 
open squares and star symbols. The results from a two-subevent method used in evaluating 
the systematic uncertainties are shown as filled circles, though they overlap the other points in 
all but the most peripheral bin. 



2.5.2 Cumulant method 

The cumulant method measures flow utilizing a cumulant expansion of multiparticle azimuthal 
correlations, without determining the orientation of the event plane. The idea is that if the par- 
ticles are correlated with the event-plane orientation, then there also exist correlations between 
them. In our analysis, we utilize two- and four-particle correlations. To calculate the cumulants 
of these correlations, from which the flow coefficient is extracted, we use a generating function 
of the multiparticle correlations in a complex plane ||42|. First, we evaluate the "integrated", or 
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reference, flow by constructir\g the corresponding generating function including afl particles 
from a broad [pj, rj) window, and averaging over the events in a given centrality class. The 
reference flow may not be corrected for tracking efficiency and should not be equated with 
the fully corrected integrated flow of the events. Then, the differential flow, i.e., the flow in 
a narrower phase-space window, either in pj or rj, is measured with respect to the reference 
flow. In the cumulant and Lee-Yang zeros methods, the reference flow serves the same pur- 
pose as the determination of the event-plane angle and the resolution correction factors in the 
event-plane method. In our analysis of V2{pt), the pj and // ranges for the reference flow are 
0.3 < Pj < 3GeV/c and \}]\ < 0.8, respectively. In the analysis of V2{^), the reference flow is 
obtained for the range 0.3 < px < 3GeV/c and |//| < 2.4 in order to maximize the resolution 
parameter, which improves as the charged hadron multiplicity M in the selected phase-space 
window increases. The transverse momentum restriction of px < 3GeV/c is imposed to limit 
the contributions from hadrons originating in jets, and thus reduce the nonflow correlations 
contributing to the measured elliptic anisotropy parameter. To avoid auto-correlations, the 
particles used for determining differential flow are not included in evaluating the reference 
flow. The generating function for the reference flow is calculated at several different points 
in the complex plane, and we then interpolate between these points. We use three values for 
the radius parameter, ro, and seven values for the polar angle, as described in Ref. Il42ll . The 
radius parameters are determined according to the detected charged particle multiplicity and 
the number of events analyzed in each centrality class. Each particle in the differential px or 
?/ bin is correlated to the particles used for the reference flow through a differential generating 
function. To account for the fact that the track reconstruction efficiency may vary across the 
chosen bin, we implement an efficiency correction that is applied as a track-by-track weight in 
the construction of the differential generating function. 

2.5.3 Lee- Yang zeros method 

The Lee-Yang zeros (LYZ) method 1431 144| for directly measuring the flow is based on mul- 
tiparticle correlations involving all particles in the event. It uses the asymptotic behavior of 
the cumulant expansion to relate the location of the zeros of a complex function to the mag- 
nitude of the integrated flow in the system. For a detector with a uniform detection efficiency 
for all particles in the chosen // and px window, it is thus possible to obtain the integrated 
flow in a simple one-step procedure. Since this is not the case for the CMS measurements pre- 
sented here, we replace the term "integrated flow" that is used in the literature describing the 
method 1431 144| with "reference flow". In the generating function (i.e., G^(zr), where r is the 
imaginary axis coordinate), the flow vector, constructed from all particles in the event, is pro- 
jected onto a direction that makes an arbitrary angle 9 with respect to the x axis. We use five 
different projection angles and then average the results over events from the same centrality 
class to reduce the statistical uncertainties. Subsequently, the minimum of the generating func- 
tion is found, and the differential flow is determined with respect to the reference flow. There 
are different ways to define the generating function that involve either a sum or a product of 
the individual particle contributions. An example is given in Fig.|5| showing how the minimum 
of the generating function is determined using both definitions. The results presented here are 
based on the product generating function. In our analysis of C2(Pt)/ the px and tj ranges for the 
reference flow are 0.3 < px < 12GeV/c and |?/| < 0.8, respectively. Since the Lee-Yang zeros 
method is less sensitive to jet-induced charged-particle correlations than the two-particle cu- 
mulant method, the px range of the tracks included in the determination of the reference flow 
is not restricted to low pj. In the analysis of V2{^), the pseudorapidity range for the particles 
included in the reference flow is extended to \ri\ < 2.4. 

The Lee-Yang zeros method is sensitive to multiplicity fluctuations. To evaluate the effects of 
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these fluctuations we perform a toy-model MC study. Event ensembles are generated sampling 
the multiplicity for each event from Gaussian distributions with mean and r.m.s. values com- 
parable to the ones measured in the centrality bins used in the measurement. For each particle, 
the corresponding pj, rj, and (p values are sampled from a realistic input V2{pT,rj) distribu- 
tion. These events are then analyzed using the same procedure as in the data and the resulting 
^lipT,^) values are compared to the input. We find that the Lee-Yang zeros method tends to 
underestimate the results if the r.m.s. of the multiplicity distribution is more than about 14% 
of the mean. In order to keep the systematic uncertainties below 2%, the data are analyzed in 
5%-wide centrality classes and then averaged and weighted with the charged-particle yield, 
to obtain results in wider centrality intervals needed for comparisons with other methods. Ef- 
ficiency corrections are implemented as a track-by-track weight in the differential generating 
function. 
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Figure 5: (Color online) An example of the modulus of the second harmonic Lee-Yang zero 
generating function G^(zr) as a function of the imaginary axis coordinate r for 9 = 0. Both the 
sum and product generating functions are shown, calculated from events with centrality 15- 
20%, |f/| < 0.8, and 0.3 < pj < 12GeV/c. An enlargement of |G^(z>)| around its first minimum 
is shown in the inset. 



2.6 Corrections to the anisotropy parameter V2 

In determining the C2(Pt) distributions for particles detected in the tracker, it is necessary to 
correct for the influence of misidentified (i.e., "fake") tracks on the measurements. As shown 



in Section 2.4 the fake-track contribution is particularly significant at low-px values and for 
pseudorapidities |//| > 1.6. Of particular concern is the observation that the fake tracks can 
carry a V2 signal at low pj similar to that of properly reconstructed (i.e., "real") tracks at a 
higher pj. Since the true V2 signal is very small at low pj, but increases at higher pj, the fake 
tracks may contribute a significant fraction of the measured V2 signal at low pj. 

Studies using a full MC CMS simulation of PbPb collisions based on the hydjet event genera- 
tor [|57| indicate that the component of the V2 signal due to fake tracks is relatively constant for 
Pt < 0.8GeV/c, where the fraction of fake tracks is largest. For higher pj, where the fraction 
of fake tracks is quite small, the value of V2 is consistent with the measured value from cor- 
rectly reconstructed tracks. This suggests the following simple correction scheme. Let Ndet(pT) 
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be the number of reconstructed tracks in a given pj bin, / the fraction of these tracks that are 
"fake," Ntrue the number of "true" tracks in the pj bin, and e the efficiency for reconstructing a 
true track in the bin. Then N^et — fNdet = ^Mrue- The fN^et fake tracks are characterized by a 
constant V2 value given by vtf^*^. The N^et — f^det tracks are characterized by V2^^^. Then 
the observed value vf^^ of V2 will be 



..obs 



{l-f)vr' + fv: 



,,fake 



(15) 



and so 



7,obs f^.fake 

This correction for the fake-track signal is only significant for pj values less than ?»lGeV/c. 
In this range, an empirical correction that results in values of V2 that are independent of the 
track selection requirements or fraction of fake tracks is applied using vt^^*^ = 1.3 {^2), where 
the yield-weighted average is performed over the transverse momentum range 0.3 to 3 GeV/c, 
folding in the efficiency-corrected spectra. This value for v^2^^ is also supported by MC studies 
using HYDJET. 

2.7 Systematic uncertainties 

2.7.1 Systematic uncertainties in the measurements of V2 

The systematic uncertainties in the measurements of V2 include those common to all methods, 
as well as method-specific uncertainties. They are evaluated as relative uncertainties and are 
reported as percentages relative to the measured V2 values. Since we are reporting the results 
on V2 for nonidentified charged particles, it is important to investigate the tracking efficiency as 
a function of particle species. The tracking efficiencies for charged pions, kaons, protons, and 
antiprotons are determined using a full simulation of CMS. Subsequently, the value of C2(Pt) 
for charged particles is obtained using different assumptions for the px^dependence of V2 and 
the transverse momentum spectra of each particle type, taking into account the corresponding 
reconstruction efficiencies. The results are compared to those obtained with the assumption of 
a particle-species-independent efficiency. The uncertainties in the charged particle V2 results 
are estimated to be < 0.5%, independent of the pj, tj, and centrality ranges. This uncertainty is 
listed as "Part, composition" in Tables |2]-p] 



Since the V2 value changes with centrality, an uncertainty in the centrality determination can 
lead to a shift in the V2 measurement. This uncertainty is evaluated by varying the value of the 
minimum-bias trigger efficiency to include (97 ± 3)% of the total inelastic cross section. The 
resulting uncertainty on V2 is of the order 1%, independent of the pj, rj, and centrality ranges. 
This uncertainty is listed as "Cent, determination" in Tables |2]-[9| 



The kinematic requirements used to select tracks can affect the efficiency of track finding and 
the relative fraction of fake tracks in an event. The requirements are varied from their default 
values to estimate the systematic uncertainty. For each set of requirements, corresponding 
corrections are obtained for the fake-track contribution to the V2 signal, as described in Sec- 
tion 2.6 For a given centrality and px range, the systematic uncertainty is estimated based on 
the stability of the V2 value after corrections for fake tracks, independent of the track selection 
requirements. The uncertainty is found to be directly related to the magnitude of the fake-track 
contribution, with the final results deemed unreliable when the fake-rate is higher than ?»20%. 
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For the results presented in Section |3| the systematic uncertainties from this source remain be- 
low 4% over the entire range of pj, t], and centrality, and are significantly below this value for 
Pt > 0.5GeV/c, centrality above 10%, and \r]\ < 1.6. The uncertainty in the efficiency correc- 
tions is evaluated by determining the efficiency based on the hydjet model, and by embedding 
simulated pions into PbPb events in data. Although the two resulting efficiencies do have dif- 
ferences, the uncertainty on the V2 value is small, at most 0.5%. Variations in the V2 results due 
to changing detector conditions throughout the data-taking period are studied by dividing the 
data into three subgroups and are found to be below 1% for all measurements. The combined 
uncertainties from the efficiency corrections, fake-track corrections, and variations in detector 
conditions are listed under "Corrections" in Tables |2]-0 

Additional studies of the systematic uncertainty are conducted for each method. In the event- 
plane method, flattening corrections are obtained using different procedures. The vertex de- 
pendence of the flattening parameters is examined and different subevent rj gaps are used in 
obtaining the resolution corrections. The uncertainties from these sources are found to be neg- 
ligible. The resolution corrections are measured with the three-subevent and two-subevent 
methods and are found to be consistent within the statistical uncertainties. The statistical un- 
certainties for the resolution correction factor are less than 1%, except for the most peripheral 
70-80% centrality events, where the statistical uncertainty reaches a value of 2%. We include 
the statistical uncertainties associated with the resolution-correction factors as part of the over- 
all systematic uncertainty on tracking efficiency and fake-track corrections for the event-plane 
method. 

In the cumulant method, we examine the numerical stability of the result when the radius 
parameter ro used in the interpolations of the generating function is increased or decreased by 
half of its central value. The effects of multiplicity fluctuations are studied for the Lee-Yang 
zeros and the cumulant methods by analyzing the events in finer 2.5% centrality bins, and by 
using a fixed number of particles chosen at random from each event in a given centrality class. 

The systematic uncertainties are smallest for the mid-rapidity region \ f]\ < 0.8, px > 0.5 GeV/c, 
and in the mid-central events (10-40%), and range from 2.0 to 4.5% for the different methods. 
At low Pj for the most central events, and in the forward pseudorapidity region where the 
fake-track contributions are larger, the uncertainties increase to 2.4-6.8%. Similarly, in the most 
peripheral events the uncertainties increase mostly due to multiplicity fluctuations and reach 
up to 3.2-7%, depending on the experimental method. A summary of the systematic uncer- 
tainties is presented in Tables |2j|9] 

Table 2: Systematic uncertainties in the measurement of V2{pt) for < 0.8 with the event- 
plane method for different pj and centrality ranges. 



Source 


Pt 
(GeV/c) 


Centrality 


- 10% 


10 - 70% 


70 - 80% 


Part, composition 


All 


0.5% 


0.5% 


0.5% 


Cent, determination 


All 


1.0% 


1.0% 


1.0% 


Corrections 


< 0.3 


4.0% 


2.0% 


3.0% 


0.3-0.5 


2.0% 


< 1.0% 


2.0% 


0.5-22.0 


< 1.0% 


< 1.0% 


2.0% 


Total 


< 0.3 


4.2% 


2.3% 


3.2% 


0.3-0.5 


2.3% 


1.5% 


2.3% 


0.5-22.0 


1.5% 


1.5% 


2.3% 
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Table 3: Systematic uncertainties in the measurement of V2{f]) for 0.3 < px < 3GeV/c with the 
event-plane method for different ?/ and centrality ranges. 



Source 


\v\ 


Centrality 


- 10% 


10 - 70% 


70 - 80% 


Part, composition 


All 


0.5% 


0.5% 


0.5% 


Cent, determination 


All 


1.0% 


1.0% 


1.0% 


Corrections 


0.0-1.6 


2.0% 


< 1.0% 


2.0% 


1.6-2.4 


4.0% 


2.0% 


3.0% 


Total 


0.0-1.6 


3.2% 


1.5% 


2.3% 


1.6-2.4 


4.2% 


2.3% 


3.2% 



Table 4: Systematic uncertainties in the measurement of V2{'Pt) for \ri\ < 0.8 with the two- 
particle cumulant method for different pj and centrality ranges. 



Source 


(GeV/c) 


Centrality 


- 5% 


5 - 10% 


10 - 70% 


70 - 80% 


Part, composition 


All 


0.5% 


0.5% 


0.5% 


0.5% 


Cent, determination 


All 


1.0% 


1.0% 


1.0% 


1.0% 


Multiplicity fluct. 


All 


0.5% 


1.5% 


4.0% 


4.0% 


ro parameter 


All 


0.2% 


0.2% 


0.2% 


0.1% 


Corrections 


0.3-0.5 


3.8% 


1.7% 


0.6% 


4.0% 


0.5-22.0 


2.9% 


2.1% 


0.6% 


3.0% 


Total 


0.3-0.5 


4.0% 


2.6% 


4.2% 


5.8% 


0.5-22.0 


3.2% 


2.9% 


4.2% 


5.1% 



2.7.2 Systematic uncertainties in tlie measurements of thie transverse momentum spec- 
tra and thie mean transverse momentum (pj) 

Several sources of systematic uncertainty are considered in obtaining the inclusive charged- 
particle transverse momentum distributions and their mean values, (pi)- These include the 
uncertainties in the efficiency and fake-track correction factors, the particle-species-dependent 
efficiency, and the uncertainty in the minimum-bias trigger efficiency. The effect on the overall 
normalization of the spectra is considered separately from the smaller effect on the shape of 
the spectra as a function of pj. To obtain the mean transverse momentum, different functional 
forms are used to extrapolate the spectra down to pj = 0, and the pj range over which the 
spectra are fitted is varied. The combined point-to-point systematic uncertainties on the spectra 



are presented as a function of pseudorapidity in Table 10 The total normalization uncertainty 
of the charged-particle spectra measured in each centrality interval is given in Table 11 The 
systematic uncertainties in the measurement of (pj) in different pseudorapidity intervals are 
summarized in Table [l2l 



3 Results 

The main results of the analysis using the four methods described above are as follows: 

• ^lipi) at mid-rapidity |//| < 0.8. 

• Integrated V2 at mid-rapidity |f/| < 0.8 and 0.3 < pj < 3GeV/c. 

• V2{r]) for 0.3 < pi < 3GeV/c. 
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Table 5: Systematic uncertainties in the measurement of V2{rj) for the range 0.3 < pj < 3GeV/c 
with the two-particle cumulant method for different r] and centrality ranges. 



Source 


\v\ 


Centrality 


5 - 10% 


10 - 70% 


Part, composition 


All 


0.5% 


0.5% 


Cent, determination 


All 


1.0% 


1.0% 


Multiplicity fluct. 


All 


1.5% 


4.0% 


ro parameter 


All 


0.2% 


0.2% 


Corrections 


0.0-1.6 


0.8% 


1% 


1.6-2.4 


1.5% 


1.6% 


Total 


0.0-1.6 


2.0% 


4.2% 


1.6-2.4 


2.4% 


4.4% 



Table 6: Systematic uncertainties in the measurement of 1^2 (Pt) for \t]\ < 0.8 with the four- 
particle cumulant method for different px and centrality ranges. 



Source 


Pt 
(GeV/c) 


Centrality 


5 - 10% 


10 - 40% 


40 - 60% 


60 - 70% 


Part, composition 


All 


0.5% 


0.5% 


0.5% 


0.5% 


Cent, determination 


All 


1.0% 


1.0% 


1.0% 


1.0% 


Multiplicity fluct. 


All 


5.0% 


3.0% 


5.0% 


5.0% 


ro parameter 


All 


2.0% 


3.0% 


1.0% 


0.1% 


Corrections 


0.3-0.5 


4.1% 


1.4% 


3.0% 


4.5% 


0.5-22.0 


2.2% 


1.2% 


1.1% 


3.0% 


Total 


0.3-0.5 


6.8% 


4.6% 


6.0% 


6.8% 


0.5-22.0 


5.9% 


4.5% 


5.3% 


6.0% 



We also measure the charged-particle transverse momentum spectra and their mean pj for the 
centrality and pseudorapidity ranges in which the flow is studied. 

The flow studies are performed in the 12 centrality classes listed in Table [l] Using these results, 
we examine the scaling of the integrated V2 with the participant eccentricity, as well as perform 
comparisons to measurements from other experiments. Centrality classes are regrouped to per- 
form these comparisons, i.e., the results of V2{pt), '02{ri), or integrated V2 obtained in the finer 
bins of centrality are averaged over wider bins, weighted using the corresponding d^N / dpj dtj 
spectra. The evolution of the measured elliptic anisotropy as a function of centrality, center-of- 
mass energy, and transverse particle density is studied. The scaling of V2{tj) in the longitudinal 
dimension is also examined through comparisons to RHIC data. 

3.1 Transverse momentum dependence of V2 

In Figs. |6]-|9| we present the measurement of V2 for charged particles as a function of transverse 
momentum at mid-rapidity, obtained by each analysis method. We use the notation z;2{EP} 
to refer to the measurement of V2 using the event-plane method, and V2{2}, V2{4i}, V2{\yZ} 
to refer to those using the two-particle cumulant, four-particle cumulant, and Lee-Yang zeros 
methods, respectively. Several trends can be observed and related to the physics processes 
dominating hadron production in different pj ranges. The value of V2 increases from central 
to peripheral collisions up to 40% centrality, as expected if the anisotropy is driven by the 
spatial anisotropy in the initial state [|T2l[l5ll2T| . The transverse momentum dependence shows 



3.1 Transverse momentum dependence of V2 



19 



Table 7: Systematic uncertainties in the measurement of V2{ri) for the range 0.3 < pj < 3GeV/c 
with the four-particle cumulant method for different rj and centrality ranges. 



Source 


\v\ 


Centrality 


5 - 10% 


10 - 40% 


40 - 70% 


Part, composition 


All 


0.5% 


0.5% 


0.5% 


Cent, determination 


All 


1.0% 


1.0% 


1.0% 


Multiplicity fluct. 


All 


5.0% 


3.0% 


5.0% 


yq parameter 


All 


2.0% 


3.0% 


1.0% 


Corrections 


0.0-1.6 


1.5% 


1.5% 


1.5% 


1.6-2.4 


1.8% 


2.1% 


1.9% 


Total 


0.0-1.6 


5.8% 


4.8% 


5.4% 


1.6-2.4 


5.8% 


5.0% 


5.6% 



Table 8: Systematic uncertainties in the measurement of ^^2(^1) for \rj\ < 0.8 with the Lee-Yang 
zeros method for different px arid centrality ranges. 



Source 


Pt 
(GeV/c) 


Centrality 


5 - 10% 


10 - 40% 


40 - 50% 


Part, composition 


All 


0.5% 


0.5% 


0.5% 


Cent, determination 


All 


1.0% 


1.0% 


1.0% 


Multiplicity fluct. 


All 


0.1% 


0.9% 


1.9% 


Corrections 


0.3-0.5 


2.5% 


1.7% 


0.7% 


0.5-22.0 


1.5% 


1.0% 


0.6% 


Total 


0.3-0.5 


2.7% 


2.2% 


2.3% 


0.5-22.0 


1.9% 


1.8% 


2.3% 



a rise of V2 up to pj ~3GeV/c and then a decrease. As a function of centrality, a tendency 
for the peak position of the C2(Pt) distribution to move to higher px in more central collisions 
is observed, with the exception of the results from the most peripheral collisions in the two- 
particle cumulant and the event-plane methods. 

In ideal hydrodynamics the azimuthal anisotropy continuously increases with increasing px l|T2l 
[21 J. The deviation of the theory from the RHIC data at px 2-3GeV/c has been attributed to 
incomplete thermalization of the high-px hadrons, and the effects of viscosity. Indeed, viscous 
hydrodynamic calculations [151 [lH |65l show that the shear viscosity has the effect of reducing 
the anisotropy at high px. At px ^ 8GeV/c, where hadron production is dominated by jet 
fragmentation, the collective-flow effects are expected to disappear [ 15ll2T| . Instead, an asym- 
metry in the azimuthal distribution of hadron emission with respect to the reaction plane could 
be generated by path-length-dependent parton energy loss l|28ll30ll3Tl . For events with simi- 
lar charged-particle suppression |66|, but different reaction-zone eccentricity, one might expect 
that the geometric information would be imprinted in the elliptic anisotropy signal. The upper 
panels of Figs. |6]-[9] (centrality 0-35%) show a trend that is consistent with this expectation. In 
more central events, where the eccentricity is smaller, the elliptic anisotropy value is system- 
atically lower. In more peripheral collisions (centrality 35-80%), there is a complex interplay 
between the reduced energy loss and the increase in eccentricity that influence the V2{pj) value 
in opposite directions. The data presented here provide the basis for future detailed compar- 
isons to theoretical models. 

An important consideration in interpreting the C2(Pt) results is the contribution from non- 
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Table 9: Systematic uncertainties in the measurement of V2{^) for 0.3 < pj < 3GeV/c with the 
Lee-Yang zeros method for different rj and centrality ranges. 



Source 


\v\ 


Centrality 


5 - 10% 


10 - 40% 


40 - 50% 


Part, composition 


All 


0.5% 


0.5% 


0.5% 


Cent, determination 


All 


1.0% 


1.0% 


1.0% 


Multiplicity fluct. 


All 


0.1% 


0.9% 


1.9% 


Corrections 


0.0-1.6 


1.3% 


1.0% 


0.8% 


1.6-2.4 


1.5% 


1.4% 


1.3% 


Total 


0.0-1.6 


1.7% 


1.8% 


2.4% 


1.6-2.4 


1.9% 


2.0% 


2.5% 



Table 10: Point-to-point systematic uncertainties in the measurement of the charged-particle 
spectra in different pseudorapidity intervals. 



Source of Uncertainty 


\ri\ < 0.8 


0.8 < \ri\ < 1.2 


2.0 < \r]\ < 2.4 


Tracking Efficiency 


5% 


8% 


13% 


Particle Composition 


1% 


1% 


2% 


Trigger Efficiency 


3% 


3% 


3% 


Total 


6% 


9% 


14% 



flow correlations and initial-state eccentricity fluctuations. To aid in assessing the magnitude 
of these effects and their evolution with the centrality of the collisions, the results of V2{pj) 



obtained by all methods at mid-rapidity are compared in Fig. 10 for 12 centrality classes. 
The four methods show differences as expected due to their sensitivities to nonflow contri- 
butions llZlllllllll and eccentricity fluctuations Il34ll35ll38l. 

The method that is most affected by nonflow correlations is the two-particle cumulant, because 
of the fact that the reference and the differential flow signals are determined in the same pseu- 
dorapidity range. The event-plane method is expected to be similarly affected if dedicated 
selections are not applied to reduce these contributions. In our analysis, the particles used in 
the event-plane determination and the particles used to measure the flow are at least 3 units of 
pseudorapidity apart, which suppresses most nonflow correlations. The differences between 
the two-particle cumulant and the event-plane methods are most pronounced at high pj and 
in peripheral collisions, where jet-induced correlations dominate over the collective flow. 

In a collision where M particles are produced, direct A:-particle correlations are typically of order 
1/ M^^^, so that they become smaller as k increases. Therefore, the fourth-order cumulant and 
the Lee-Yang zeros methods are expected to be much less affected by nonflow contributions 
than the second-order cumulant method Il37il42ll43l| . This trend is seen in our data. 

3.2 Centrality dependence of integrated V2 and eccentricity scaling 

To obtain the integrated V2 values as a function of centrality at mid-rapidity, ^2(^1) measure- 
ments are averaged over pj, weighted by the corresponding charged-particle spectrum. The 
integration range 0.3 < pj < 3GeV/c is limited to low pj to maximize the contribution from 
soft processes, which facilitates comparisons to hydrodynamic calculations. 



The centrality dependence of the integrated Vz at mid-rapidity |/y | < 0.8 is presented in Fig. 11 



for the four methods. The V2 values increase from central to peripheral collisions, reaching 



3.2 Centrality dependence of integrated C2 and eccentricity scaling 
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Table 11: Normalization uncertainty in the measurement of the charged-particle spectra in dif- 
ferent centrality intervals resulting from the uncertainty in the minimum-bias trigger efficiency. 



Centrality Range 


Normalization Uncertainty 


0-5% 


0.4% 


5-10% 


1.0% 


10-15% 


1.7% 


15-20% 


2.3% 


20-25% 


3.1% 


25-30% 


4.1% 


30-35% 


5.0% 


35-40% 


6.1% 


40-50% 


8.0% 


50-60% 


12% 


60-70% 


16% 


70-80% 


21% 



Table 12: Systematic uncertainty of the mean pj of charged particles from each source and in 
total as a function of pseudorapidity. 



Source of Uncertainty 


\f]\ < 0.4 


0.8 < \r]\ < 1.2 


2.0 < < 2.4 


Fit function 


3% 


3% 


4% 


Trigger Efficiency 


1.5% 


1.5% 


1.5% 


Tracking Efficiency 


2% 


2% 


2.5% 


Total 


3.9% 


3.9% 


4.9% 



a maximum in the 40-50% centrality range. In the more peripheral collisions, a decrease in 
V2 is observed in the event-plane and four-particle cumulant measurements, while the values 
obtained with the two-particle cumulant method remain constant within their uncertainties. 
The results for V2{2} are larger than those for C2{EP}, while the f 2{4} and f 2{LYZ} values are 
smaller. To facilitate a quantitative comparison between the methods, including their respec- 



tive systematic uncertainties, the bottom panel of Fig. 11 shows the results from the cumulant 
and the Lee-Yang zeros methods divided by those obtained from the event-plane method. The 
boxes represent the systematic uncertainty in the ratios, excluding sources of uncertainty com- 
mon to all methods. The ratios are relatively constant in the 10-60% centrality range, but the 
differences between the methods increase for the most central and the most peripheral colli- 
sions. These findings are similar to results obtained by the STAR experiment at RHIC ||67|. 
Below, we further investigate the differences in the V2 values returned by each method. 

The collective motion of the system, and therefore the anisotropy parameter, depend on the 
initial shape of the nucleus-nucleus collision area and the fluctuations in the positions of the 
interacting nucleons. By dividing V2 by the participant eccentricity, one may potentially remove 
this dependence across centralities, colliding species, and center-of-mass energies, enabling a 
comparison of results in terms of the underlying physics driving the flow. 



In Fig. 12 we examine the centrality dependence of the eccentricity-scaled anisotropy parame- 
ter obtained with the event-plane and cumulant methods at mid-rapidity < 0.8. The partic- 
ipant eccentricity and its cumulant moments are obtained from a Glauber-model simulation, as 



discussed in Section 2.3 The statistical and systematic uncertainties in the integrated V2 mea- 
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Figure 6: (Color online) Results from the event-plane (EP) method for V2 as a function of pj at 
mid-rapidity |?/ 1 < 0.8 for the 12 centrality classes given in the legend. The error bars show the 
statistical uncertainties only. 



surements are added in quadrature and represented by the error bars. The dashed lines show 
the systematic uncertainties in the eccentricity determination. In the left panel of Fig. 12 the re- 
sults from each method are divided by the participant eccentricity, epart- The data show a near- 
linear decrease from central to peripheral collisions, with differences between methods that 
were already observed in Fig. 11 In the right panel, the V2 values for the cumulant measure- 
ments are scaled by their respective moments of the participant eccentricity, thus taking into 
account the corresponding eccentricity fluctuations [34-38 1. With this scaling, the two-particle 
cumulant and the event-plane results become nearly identical, except for the most central and 
the most peripheral collisions, where the cumulant results are more affected by nonflow con- 
tributions. This is expected [37J because in our application of the method there is no separation 
in rapidity between the particles used for the reference flow and those used in the differential 
flow measurement. In the centrality range of 15-40%, the four-particle cumulant measurement 
of C2{4}/e{4} is also in better agreement with the other two methods. This indicates that the 
main difference in the results from the different methods could be attributed to their sensitiv- 
ity to eccentricity fluctuations. In the most central events, where the eccentricity e{4} is very 
small, and in the most peripheral events, where the fluctuations are large, t?2{4}/e2{4} devi- 
ates from the common scaling behavior. For centralities above 50% the differences between the 
four-particle cumulant method and the other two methods do not seem to be accounted for 
by the initial-state fluctuations, as described in our implementation of the Glauber model. In 
this centrality range, the event-by-event fluctuations in the eccentricity are non-Gaussian due 
to the underlying Poisson distributions from discrete nucleons (3411381 and are more difficult to 
model. It has also been suggested [,34, ,37] that when the event-plane resolution is smaller than 
~0.6, as is the case for the peripheral collisions studied in CMS, the results from the event-plane 
method should be evaluated using the two-particle cumulant eccentricity e{2}, rather than the 
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Figure 7: (Color online) Second-order cumulant results for as a function of px at mid-rapidity 
I?/ 1 < 0.8 for the 12 centrality classes given in the legend. The error bars show the statistical 
uncertainties only. 



participant eccentricity epart- We have used a common definition of eccentricity (Cpart) for all 
centrality classes studied in our event-plane analysis. This would lower the measurements of 
z;2{EP}/e by about 10% in the most peripheral collisions, which is not sufficient to reconcile 
the differences between the event-plane and the four-particle cumulant results. 

Another model of the initial state that has been used in the literature l|38ll39ll68ti7T| . but has not 
been explored here, is the color glass condensate (CGC) model |I72I| , which takes into account 
that at very high energies or small values of Bjorken x, the gluon density becomes very large 
and saturates. The CGC model predicts eccentricities that exceed the Glauber-model eccentric- 
ities by an approximately constant factor of around 1.2, with some deviation from this behavior 
in the most central and most peripheral collisions ||38l[70| . The results presented here may give 
further insight into the nature of the initial-state fluctuations, especially in the regions where 
the eccentricity fluctuations become non-Gaussian. 

3.3 Pseudorapidity dependence of 

The pseudorapidity dependence of the anisotropy parameter provides additional constraints 
on the system evolution in the longitudinal direction. To obtain the V2{y]) distribution with 
the event-plane method, we first measure v^iji-r) in pseudorapidity bins of A// = 0.4, and then 
average the results over the range 0.3 < pj < 3GeV/c, weighting with the efficiency and fake- 
rate-corrected spectrum. 

For the cumulant and Lee-Yang zeros methods, the measurements are done using all particles 
in the range \r]\ < 2.4, and either 0.3 < px < 3GeV/c or 0.3 < pj < 12GeV/c in the generating 
function, to obtain the reference flow, and then extracting the pseudorapidity dependence in 
small pseudorapidity intervals of Ar] = 0.4. Tracking efficiency and fake-rate corrections are 
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Figure 8: (Color online) Fourth-order cumulant results for V2 as a function of px at mid-rapidity 
< 0.8 for the 10 centrality classes given in the legend. The error bars show the statistical 
uncertainties only. 
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Figure 9: (Color online) Lee-Yang zeros results for V2 as a function of pj at mid-rapidity 
I?/ 1 < 0.8 for the 8 centrality classes given in the legend. The error bars show the statistical 
uncertainties only. 
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Figure 10: (Color online) Comparison of the four different methods for determining V2 as a 
function of pj at mid-rapidity (|f/| < 0.8) for the 12 centrality classes given in the figures. The 
error bars show the statistical uncertainties only. 



applied using a track-by-track weight in forming the differential generating functions. As a 
crosscheck, we have confirmed that at mid-rapidity the values obtained with this method agree 
with the ones obtained from a direct yield-weighted average of the 1^2(^1) results from Figs.j^ 
|9| within the stated systematic uncertainties. 



As observed at mid-rapidity ( I ;y I < 0.8)inFig.ll the values of t?2{4} and ?72{LYZ} are in agree 



ment and are smaller than C2{2} and 172{EP}. This behavior persists at larger pseudorapidity. 



as shown in Fig. 13 which suggests that similar nonflow correlations and eccentricity fluctua- 
tions affect the results over the full measured pseudorapidity range. The results show that the 
value of V2{ti) is greatest at mid-rapidity and is constant or decreases very slowly at larger val- 
ues of \tj\. This behavior is most pronounced in peripheral collisions and for the two-particle 
cumulant method, which is most affected by nonflow contributions. 

To assess whether the observed decrease in V2{Tj) in the forward pseudorapidity region in pe- 
ripheral PbPb collisions is due to a pseudorapidity dependence in the 1^2 (Pt) distributions or 



in the underlying charged-particle spectra, in Fig. 14 we examine the values of 1^2 (Pt) obtained 



with the event-plane method for several pseudorapidity intervals in each of the 12 central- 



ity classes shown in Fig. 13 From the most central events up to 35-40% centrality there is 
no change in the 1^2 (Pt) distributions with pseudorapidity within the statistical uncertainties. 
Therefore, any change in the V2{yi) distribution can be attributed to changes in the underlying 
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Figure 11: (Color online) Top panel: Integrated V2 as a function of centrality at mid-rapidity 
< 0.8 for the four methods. The boxes represent the systematic uncertainties. The mag- 
nitudes of the statistical uncertainties are smaller than the size of the symbols. Bottom panel: 
The values from three of the methods are divided by the results from the event-plane method. 
The boxes represent the systematic uncertainties excluding the sources that are common to all 
methods. The magnitudes of the statistical uncertainties are smaller than the size of the sym- 
bols. 



charged-particle transverse momentum spectra. A gradual decrease is observed in the C2(Pt) 
values at forward pseudorapidity (2.0 < < 2.4) in more peripheral events. For the 70- 
80% centrality class the values of 1^2(^1) decrease by approximately 10% between the central 
pseudorapidity region \ri\ < 0.4 and the forward region 2.0 < \rj\ < 2.4. Thus, the pseudora- 



pidity dependence in V2{f]) for peripheral collisions observed in Fig. 13 is caused by changes in 



the V2{pt) distributions with pseudorapidity, as well as changes in the underlying transverse 



momentum spectra presented in Section 3.4 



3.4 Centrality and pseudorapidity dependence of the transverse momentum 
distributions 

Elliptic flow measures the azimuthal anisotropy in the invariant yield of the final-state par- 
ticles. Therefore, the charged-particle transverse momentum distributions influence the ob- 
served results. The soft-particle-production mechanism and the evolution of the expanding 
nuclear medium are reflected in the low-px range of the transverse momentum spectra. In the 
hydrodynamics calculations, measurements of the pseudorapidity density of charged particles 
produced in collisions with different centrality constrain the description of the initial entropy 
and the energy density distribution in the collision zone, while the mean transverse momentum 
of the particle spectra constrains the final temperature and the radial-flow velocity of the sys- 
tem. With additional input on the equation-of-state that is typically provided by lattice QCD, 
the hydrodynamics calculations then provide a description of the system evolution from some 
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Figure 12: (Color online) Left panel: Centrality dependence of the integrated V2 divided by the 
participant eccentricity, epart, obtained at mid-rapidity < 0.8 for the event-plane and two- 
and four-particle cumulant methods. Right panel: The measurements of V2 as a function of 
centrality are the same as in the left panel, but here the two-particle and four-particle cumulant 
results are divided by their corresponding moments of the participant eccentricity, e{2} and 
e{4}. In both panels, the error bars show the sum in quadrature of the statistical and systematic 
uncertainties in the V2 measurement, and the lines represent the systematic uncertainties in the 
eccentricity determination. 



initial time, when local thermal equilibrium is achieved, to the "freeze-out", when the particle 
interactions cease. 

We have measured the charged-particle transverse momentum spectra for 12 centrality classes 
over the pseudorapidity range |?/| < 2.4 in bins of At] = 0.4. Examples of these distributions 
for the mid-rapidity < 0.4) and forward-rapidity (2.0 < |?/| < 2.4) regions are shown in 



Fig. 15 These results extend the measurements of charged-particle spectra previously reported 
by CMS [66 J down to pj = 0.3GeV/c and to forward pseudorapidity. The measurements pre- 
sented here are in good agreement with the results in Ref. II66I in their common ranges of pj, 
pseudorapidity, and centrality. 

The evolution of the charged-particle spectra with centrality and pseudorapidity can be quan- 
tified in terms of the mean pj of the transverse momentum distributions. The values of (pj) 



as a function of Npart/ which is derived from the centrality of the event, are shown in Fig. 16 
In each pseudorapidity interval, the values of {pj) increase with Npart up to Npart ~150 anc 
then saturate, indicating that the freeze-out conditions of the produced system are similar over 
a broad range of collision centralities (0-35%). This behavior is in contrast to the centrality 



dependence in the integrated V2 values at mid-rapidity shown in Fig. 11 that vary strongly in 
this centrality range. On the other hand, in more peripheral collisions (centrality greater than 
35%) the V2{pt) values do not vary much with centrality at low pj, as shown in the bottom 
panels of Figs. |6}|9| while the spectral shapes change, as indicated by the (px) measurement. 



This behavior is qualitatively similar at forward pseudorapidities, as shown in Figs. 14 and 16 
These measurements taken together will help in understanding the early-time dynamics in the 
system evolution, which is reflected in the elliptic flow, and the overall evolution through the 
hadronic stage, which is reflected in the charged-particle spectra. 
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Figure 13: (Color online) Pseudorapidity dependence of V2 for 0.3 < px < 3 GeV/c with all four 
methods in 12 centrality classes. The boxes give the systematic uncertainties. The magnitudes 
of the statistical uncertainties are smaller than the size of the symbols. 



3.5 Comparison with other measurements of V2 at the LHC 

Results on the elliptic anisotropy measured in PbPb collisions in a/snn = 2.76 TeV have previ- 
ously been reported by the ALICE [13J and ATLAS [14J experiments. A comparison of V2{2} 
and V2{4:} as a function of pj in the 40-50% centrality class for |?/| < 0.8 from CMS and AL- 



ICE (T3\ is shown in Fig. 17 The error bars give the statistical uncertainties, and the boxes 
represent the systematic uncertainties in the CMS measurements. The two measurements are 
in good agreement over their common pj range. 

A comparison of 1^2 (Pt) obtained with the event-plane method at mid-rapidity from CMS and 



ATLAS ||T4| is presented in Fig. 18 for the centrality ranges of the ATLAS measurement. The 
error bars show the statistical and systematic uncertainties added in quadrature.The results are 
in good agreement within the statistical and systematic uncertainties. 



4 Discussion 

We compare the CMS elliptic flow measurements presented in Section |3] with results obtained 
at RHIC by the PHENIX, STAR, and PHOBOS experiments. Since each method for measur- 
ing V2 has a different sensitivity to nonflow correlations and initial-state fluctuations, we make 
these comparisons for measurements conducted with the same method and with similar kine- 



matic requirements in the method's implementation. In Fig. 19 the mid-rapidity measurement 
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Figure 14: (Color or\line) Results from the everit-plane method for 1^2 (Pt) in three pseudora- 
pidity regions and in 12 centrality classes. The error bars show the statistical uncertainties that 
in most cases have magnitudes smaller than the size of the symbol. 



of V2{pj) with the event-plane method in CMS is compared to results from PHENIX ||731 for 
\/snn = 200 GeV AuAu collisions. For the PHENIX measurement, the event plane was de- 
termined at forward pseudorapidities, \rj\ = 3.1-3.9, while the 1^2 (Pt) measurement was per- 
formed in the pseudorapidity interval \ri\ < 0.35, thus providing a separation of at least 2.75 
units of pseudorapidity between the charged particles used for the C2(Pt) analysis and the 
particles used in the event-plane determination. This procedure is comparable to the CMS 
approach, where a separation of at least 3 units of pseudorapidity is used. These large pseudo- 
rapidity gaps are expected |37| to suppress nonflow contributions in both measurements. The 
pseudorapidity interval for the CMS measurement is wider, \rj\ < 0.8, but since the pseudo- 
rapidity dependence of V2{j]) was shown to be weak (see Fig. 13 , this difference should not 
influence the comparison of the results. The top panels in Fig. 19 show the measurements of 
172 (Pt) from CMS (closed symbols) and PHENIX (open symbols) for several centrality classes. 
The error bars represent the statistical uncertainties only. The shape of the 02(Pt) distributions 
and the magnitude of the signals are similar, in spite of the factor of ~14 increase in the center- 
of-mass energy. To facilitate a quantitative comparison of these results, the CMS measurements 
are fitted with a combination of a fifth-order polynomial function (for pj < 3.2GeV/c) and a 
Landau distribution (for 3 < px < 7GeV/c). There is no physical significance attributed to 
these functional choices, other than an attempt to analytically describe the CMS C2(Pt) distri- 
butions, so that the value of V2 can be easily compared to results from other experiments, which 
have been obtained with different px binning. The results from the fits are plotted as solid lines 
in the top panels of Fig. 19 In the bottom panels, the fit function is used to evaluate 1^2 (px) 



30 



4 Discussion 



f 

10^ 

1 



o 10 



Q. 

T3 



10 
T3 

^ 10"^ 



MO"" 
10-^° 



-2 _ 



10 



12 



CMS PbPb 



I I I I I I I I I 

0-5% (x1 0^) 
5-10% (x10\ 



\^ = 2.76TeV ^ io-i5%(x10^) j[ 



1 5-20% (x 10^) 



|2\ ~- 




I I I I I I I I I I I I I I I I I I I 



2.0<|ri|<2.4 



I I I I 

A 30-35% (x 10") ■ 
« 35-40% (xioi 

* 40-50% (x 10"^) . 

* 50-60% (x 10"^) . 

5. 




I I I I I I I I I I I I I I I I I I I 



I I I I 







2 3 
(GeV/c) 







2 3 
(GeV/c) 



Figure 15: (Color online) Inclusive charged-particle spectra at mid-rapidity (left) and forward 
rapidity (right), for the 12 centrality classes given in the legend. The distributions are offset by 
arbitrary factors given in the legend for clarity. The shaded bands represent the statistical and 
systematic uncertainties added in quadrature, including the overall normalization uncertain- 
ties from the trigger efficiency estimation. 



at the Pt values for each data set, and then to form the ratios between the CMS fit values and 
the PHENIX data, and the CMS fit to the CMS measurements. The error bars represent the sta- 
tistical uncertainties. The systematic uncertainties from the CMS and PHENIX measurements 
are added in quadrature and plotted as shaded boxes. The C2(Pt) values measured by CMS 
are systematically higher than those from PHENIX in all centrality classes and over the entire 
transverse momentum range measured by PHENIX. The relative deviations of are the order 
10%, except for the most peripheral collisions where they reach 15%. 

A similar comparison is carried out for the two-particle and four-particle cumulant methods. 



In Fig. 20 the results from the STAR experiment [74J for AuAu collisions at i/s]vn = 200 GeV 
in the 20-60% centrality range are compared with the CMS measurements. The pseudorapid- 
ity interval for the STAR measurement is \rj\ < 1.3, compared to |//| < 0.8 for CMS. These t] 
ranges are within a pseudorapidity region in which the V2{ri) values only weakly depend on 
the pseudorapidity. The kinematic selections imposed on the charged particles used in deter- 
mining the reference flow are also similar in the CMS and STAR measurements. The top panels 



of Fig. 20 show the filpx) distributions for the two-particle (left) and four-particle (right) cu- 



mulant method from both experiments, along with fits to the CMS data (lines). The functional 



form used for the fit of the CMS i?2(Pt) distributions is the same as the one used in Fig. 19 The 
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Figure 16: Mean transverse momentum of the charged-particle spectra as a function of A/part 
in three pseudorapidity intervals marked in the figure. The error bars represent the quadratic 
sum of the statistical and systematic uncertainties. 



bottom panels in Fig. 20 show the ratios of the fits to the CMS data to the actual measurements 
from CMS and STAR. The error bars represent the statistical uncertainties. The systematic un- 
certainties from the CMS and STAR measurements are added in quadrature and plotted as 
shaded boxes. At low pj, the C2 (Pt) values measured by CMS are larger than in the STAR data, 
but the relative deviations are smaller than 5% for the four-particle cumulant method, and are 
of the order 10-15% for the two-particle cumulant method. Taken together, the comparisons to 
the RHIC results in Figs. 19 and 20 indicate only a moderate increase in t72(PT) at low pj from 
the highest RHIC energy to the LHC, despite the large increase in the center-of-mass energy. 



In Fig. 21 we examine the ^/s^n dependence of the integrated V2 from mid-central collisions 
spanning y^Snn = 4.7 GeV to ^s^n = 2.76 TeV. The CMS measurement is obtained with the 
event-plane method in the 20-30% centrality class by extrapolating the viipj) and the charged- 
particle spectra down to pj = 0. In the extrapolation it is assumed that V2{0) = 0, and the 
charged-particle yield is constrained to match the dNch/ drj values measured by CMS [75 J. The 
low-energy data are from Refs. 1201 HSl SH I76ti8l1, as compiled in Ref. [79J and tabulated in 
Ref. [76] . The error bars for the low-energy data represent the statistical uncertainties. For the 
CMS data the error bar is the quadrature sum of the statistical and systematic uncertainties. 
The integrated V2 values increase approximately logarithmically with a/sjvn over the full en- 
ergy range, with a 20-30% increase from the highest RHIC energy to that of the LHC. This 
has contributions from the increase in the mean pj of the charged-particle spectra with -^s^jv, 
shown in Fig. 22 and from the moderate increase in the V2{pt) distributions at low pj, shown in 



Fig. 19 We note that the centrality selections, the collision species, and the methods employed 
in the integrated V2 measurements are not identical in all experiments, so the comparison pre- 



sented in Fig. 21 is only approximate. Further comparisons to results from lower energies are 



presented in Figs. 23 - 25 



In ideal hydrod5mamics, the eccentricity-scaled elliptic flow is constant over a broad range 
of impact parameters; deviations from this behavior are expected in peripheral collisions, in 
which the system freezes out before the elliptic flow fully builds up and saturates [32 J. A 
weak centrality and beam-energy dependence is expected through variations in the equation- 
of-state. In addition, the system is also affected by viscosity, both in the sQGP and hadronic 
stages 11221 |68[|83ll84| of its evolution. Therefore, the centrality and a/snn dependence ofv2/e 
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Figure 17: (Color online) The values of C2{2} and 1^2 {4} obtained with the cumulant method, as 
a function of pj from CMS (closed symbols) and ALICE (open symbols) ||T3| . measured in the 
range \ri\ < 0.8 for the 40-50% centrality class. The error bars show the statistical uncertainties, 
and the boxes give the systematic uncertainties in the CMS measurement. 
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Figure 18: (Color online) Comparison of results for 02(Pt) obtained with the event-plane 
method from CMS (closed symbols) and ATLAS (open symbols) for the centrality classes 
marked in the figure. The error bars show the statistical and systematic uncertainties added 
in quadrature. 
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Figure 19: (Color online) Top panels: Comparison of 1^2(^1) using the event-plane method 
as measured by CMS (solid circles) at ^ysMN = 2.76 TeV, and PHENIX [73] (open diamonds) 
at a/S]vn = 200 GeV for mid-rapidity {\r]\ < 0.8 and \r]\ < 0.35, respectively). The error bars 
represent the statistical uncertainties. The solid line is a fit to the CMS data. Bottom panels: 
Ratios of the CMS fit to the PHENIX data (open diamonds) and to the CMS data (open circles). 
The error bars show the statistical uncertainties, while the shaded boxes give the quadrature 
sum of the CMS and PHENIX systematic uncertainties. 



can be used to extract the ratio of the shear viscosity to the entropy density of the system. 



In Fig. 23 the integrated V2 obtained from the event-plane method is divided by the eccentric- 
ity of the collisions and plotted as a function of Npart/ which is derived from the centrality of 
the event. The result is compared to lower-energy AuAu and CuCu measurements from the 
PHOBOS experiment [,34 1. For the CMS measurement, the value of V2 is divided by the partici- 
pant eccentricity Cpart since the event-plane resolution factor shown in Fig. His greater than 0.6 
for all but the most central and most peripheral event selections in our analysis. It has been ar- 
gued l|34[|37l that for lower-resolution parameters, the event-plane method measures the r.m.s. 
of the azimuthal anisotropy, rather than the mean, and therefore, the relevant eccentricity pa- 



rameter in this case should be the second-order cumulant eccentricity e{2} = \J {e^^^^). Thus, 

the comparison with the PHOBOS V2 results, which were obtained with low event-plane reso- 
lution, is done by implementing this scaling using the data from Ref. Il34l . An approximately 
25% increase in the integrated V2 scaled by the eccentricity between RHIC and LHC energies is 
observed, and with a similar Npart dependence. 

It was previously observed I l34l [79l [85^1 that the v^/e values obtained in different collision sys- 
tems and varying beam energies scale with the charged-particle rapidity density per unit trans- 
verse overlap area (1/ S) (dNch/ di/), which is proportional to the initial entropy density. In ad- 
dition, it has been pointed out f69] that in this representation the sensitivity to the modeling of 
the initial conditions of the heavy-ion collisions is largely removed, thus enabling the extrac- 
tion of the shear viscosity to the entropy density ratio from the data through the comparison 
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Figure 20: (Color online) Top panels: Comparison of V2{pj) using the two-particle (left) and the 
four-particle (right) cumulant method as measured by CMS (solid circles) at ^/sj^ = 2.76 TeV, 
and STAR [74J (open stars) at ^/s^ = 200 GeV at mid-rapidity < 0.8 and \r]\ < 1.3, 
respectively). The error bars represent the statistical uncertainties. The line is a fit to the CMS 
data. Bottom panels: Ratios of the CMS fit values to the STAR data (open diamonds) and to 
the CMS data (open circles). The error bars show the statistical uncertainties, while the shaded 
boxes give the quadrature sum of the CMS and STAR systematic uncertainties. 



with viscous hydrodynamics calculations. With the factor of 2.1 increase in the charged-particle 
pseudorapidity density per participant pair, (dN^h/ d;/)/ (Npart/2), from the highest RHIC en- 
ergy to the LHC 11751 [86l , this scaling behavior can be tested over a much broader range of 
initial entropy densities. In Fig. 24 we compare the CMS results for V2/e from the event- 
plane method to results from the PHOBOS experiment |34| for CuCu and AuAu collisions 
with ^/saja/^ = 62.4 GeV and 200 GeV. 

At lower energies, the scaling has been examined using the charged-particle rapidity density 
dNch/ dy IMl IZH |85|. However, since we do not identify the species of charged particles in 
this analysis, we perform the comparison using (1/ S) [dN^h/ drj) to avoid introducing uncer- 
tainties related to assumptions about the detailed behavior of the identified particle transverse 
momentum spectra that are needed to perform this conversion. In Fig. |24j the charged-particle 
pseudorapidity density dNch/dtj measured by CMS Il75| is used, and the value of the inte- 
grated V2 for the ranges < px < 3GeV/c and < 0.8. The transverse nuclear-overlap area S 
and the participant eccentricity are listed in Table [T] The PHOBOS results from Ref . {3^ used 
dNch/ dy, and applied two factors to perform the Jacobian transformation from dN^h/ dt]: the x 
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Figure 21: (Color online) The CMS integrated V2 values for 20-30% centrality from the range 
\ri\ < 0.8 and < px < 3GeV/c obtained using the event-plane method is compared as a 
function of a/snn to results at mid-rapidity and similar centrality from ALICE |13|, STAR [76 1, 
PHENIX [23 1, PHOBOS [49, 77, 78\, NA49 LZ9[, E877 [H, and CERES [81[. The error bars for the 
lower-energy results represent statistical uncertainties; for the CMS and ALICE measurements 
the statistical and systematic uncertainties are added in quadrature. 
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Figure 22: Mean transverse momentum of the charged-particle spectra as a function of Npart 
measured by CMS in PbPb collisions at a/sjvn = 2.76 TeV (closed circles) and by STAR ||82| in 
AilAu collisions at i/snn = 200 GeV (open circles). The error bars represent the quadratic sum 
of statistical and systematic uncertainties. 
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Figure 23: (Color online) The CMS integrated V2 values from the event-plane method divided 
by the participant eccentricity as a function of Npart with \r]\ < 0.8 and < px < 3GeV/c. 
These results are compared with those from PHOBOS fMI for different nuclear species and 
collision energies. The PHOBOS V2 values are divided by the cumulant eccentricity £{2} (see 
text). The error bars give the statistical and systematic uncertainties in the V2 measurements 
added in quadrature. The dashed lines represent the systematic uncertainties in the eccentricity 
determination. 



axis was scaled by a factor 1.15, and the y axis by 0.9. Both of these factors are reversed in order 



to compare the CMS and PHOBOS measurements in Fig. 24 As in Fig. 23 the PHOBOS data 
are scaled by e{2}, while the CMS data are scaled by the participant eccentricity epart, taking 
into account the event-plane resolution factors in the two measurements. The CMS result ex- 
tends to very peripheral collisions (70-80% centrality), which allows for a significant overlap in 
the transverse charged-particle density measured at RHIC and LHC. Despite the large system- 
atic uncertainties quoted for the PHOBOS measurements, the data are in good agreement over 
the common (l/S)(dNch/ drj) range. A smooth increase in V2/ £ proportional to the transverse 
particle density is observed over the entire measured range, except for a small decrease in the 
most central collisions in both the RHIC and LHC data. The theoretical predictions [ 68117111831 
for the ^Snn dependence of the transverse-particle-density scaling of differ, and do not 
generally predict a universal behavior. The data presented here provide constraints on the 
model descriptions of the dynamical evolution of the system, and thus should aid the reliable 
extraction of the transport properties of the hot QCD medium from data. 

The V2{rj) results can be used to test theoretical descriptions of the longitudinal d}mamics in 
the expanding system, as they have been shown to be sensitive to the choice of the initial con- 
ditions, the event-by-event fluctuations in the eccentricity, and the viscosity in the sQGP and 
the hadronic stages of the system evolution I|39ll87| . The PHOBOS experiment observed 1411 
that the elliptic flow measured over a broad range of collision energies (a/smn = 19.6, 62.4, 
130, and 200 GeV) exhibited longitudinal scaling extending over several units of pseudorapid- 
ity when viewed in the rest frame of one of the incident nuclei. A similar phenomenon for 
soft-particle yields and spectra is known as limiting fragmentation II88I . Furthermore, with 
increasing i/s^N/ this beam-energy independence of V2{f]) was found to extend over an in- 
creasingly wider pseudorapidity range [41 J. To investigate the potential continuation of ex- 
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Figure 24: (Color online) Eccentricity-scaled V2. as a function of the transverse charged-particle 
density from CMS and PHOBOS ||34| . The error bars include both statistical and systematic 
uncertainties in V2. The dashed lines represent the systematic uncertainties in the eccentricity 
determination. 



tended longitudinal scaling of the elliptic flow to LHC energies, in Fig. 25 we compare the 
pseudorapidity dependence of C2('/) measured by CMS with that measured by PHOBOS at 
y'sjvN = 200 GeV in three centrality intervals. Neither the CMS nor the PHOBOS measure- 
ments are performed using identified particles. The pseudorapidity {^ ) of the particles 
in the rest frame of the nuclei moving in the positive (negative) direction is approximated by 
rj^ = T] ± 1/beam/ where rj is the pseudorapidity of the particles in the center-of-mass frame, and 
1/beam = arccosh{Ei^\,/ Am j^c-^) In (i/snn [GeV] ), with Ejab denoting the energy of the beam in 
the laboratory frame, A the nuclear mass number, and the nucleon mass. In Fig. 25 the left 
(right) half of each plot depicts V2 in the rest frame of the beam moving in the positive (negative) 
lab direction. The PHOBOS V2{ri) results are from the hit-based analysis from Fig. 4 in Ref. [|77|. 
The CMS results are obtained with the event-plane method in 0.4-unit-wide bins of pseudora- 
pidity by averaging the corresponding 1^2 (Pt) distributions over the range < pj < 3GeV/c 
and folding in the efficiency-corrected charged-particle spectra. A comparable centrality inter- 
val (2.5-15% in CMS, and 3-15% in PHOBOS) is analyzed for the most central collisions, while 
for mid-central (15-25%) and more peripheral (25-50%) collisions, the centrality selections are 
the same in both experiments. In Fig. 25 the statistical uncertainties are shown as error bars, 
while the systematic ones are represented by the shaded boxes surrounding the points. The 
CMS data cover 4.8 units of pseudorapidity, but do not overlap in pseudorapidity with the 
PHOBOS data when plotted in the rest frames of the colliding nuclei. The CMS results show 
weaker pseudorapidity dependence than observed in the PHOBOS measurement. The data 
suggest a nearly boost-invariant region that is several units wide for central events but con- 
siderably smaller for peripheral ones. It has been noted |89| that if the QCD matter produced 
at mid-rapidity at RHIC is in local equilibrium, then deviations from the triangular shape of 
V2{rj) observed in the PHOBOS measurements [41] would be expected around mid-rapidity at 
LHC energies. Detailed comparisons of theoretical calculations to the results presented here 
can give new insights into the nature of the matter produced at both RHIC and LHC energies. 



38 



5 Summary 



0.08 



0.06 



0.04 



0.02 



□ PHOBOSAuAu(s^ = 200GeV 



CMS 2.5-15% 
PHOBOS 3-15% 




ITT 




25-50% 



5 

Ti + y, 







beam 



Tl-y 



beam 



beam 



beam 



11 + y, 



beam 



Tl-y, 



beam 



Figure 25: (Color online) Measurements of V2 as a function of the pseudorapidity of particles in 
the rest frame of the colliding nuclei rj ± i/beam from CMS (closed symbols) and PHOBOS (open 
symbols) [77] in three centrality intervals. The error bars show the statistical uncertainties and 
the boxes give the systematic ones. 



5 Summary 

Detailed measurements of the charged-particle azimuthal anisotropies in a/snn = 2.76 TeV 
PbPb collisions and comparisons to lower collision energy results have been presented. The 
results cover a broad kinematic range: 0.3 < px < 20GeV/c, \ri\ < 2.4, and 12 centrality 
classes from to 80%. The measurements employ four different methods that have different 
sensitivities to fluctuations in the initial conditions and nonflow correlations. The systematic 
comparison between the methods provides the possibility to explore the underlying physics 
processes that cause these differences. 

The elliptic anisotropy parameter V2{pt) for \f]\ < 0.8 is found to increase with pj up to pj ~ 
3GeV/c, and then to decrease in the range 3 < px < lOGeV/c. For transverse momenta of 
10 < Px < 20GeV/c, no strong dependence of V2 on px is observed. The study of the high-px 
azimuthal anisotropy in charged-particle production may constrain the theoretical descriptions 
of parton energy loss and its dependence on the path-length traveled through the medium. 

The shapes of the i?2(Pt) distributions are found to be similar to those measured at RHIC. At 
low Px, only a moderate increase (5-15%) is observed in the comparison between results ob- 
tained at the highest RHIC energy and the LHC, despite the large increase in the center-of-mass 
energy. The integrated V2 at mid-rapidity and in mid-central collisions (20-30% centrality) in- 
creases approximately logarithmically with a/snn- An increase by 20-30% from the highest 
RHIC energy to that of the LHC is observed, which is mostly due to the increase in the mean 
Px of the underlying charged-particle spectra. The integrated V2 signal increases from the most 
central collisions to the 40-50% centrality range, after which a decrease is observed. Conversely, 
the values of (px) increase with Npart up to Npart ~150 (from the most peripheral collisions up 
to centrality~35%) and then saturate, indicating similar freeze-out conditions in the more cen- 
tral collisions. The different methods of measuring V2 give consistent results over a broad range 
of centrality, when scaled by their respective participant eccentricity moments. Deviations from 
this scaling are observed in the most central collisions and in peripheral (centrality above 50%) 
collisions. The eccentricity-scaled V2 at mid-rapidity is measured to be approximately linear 
in the transverse particle density, and a universal scaling is observed in the comparison of re- 
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suits from different collision systems and center-of-mass energies measured at RHIC and the 
LHC. The value of V2{f]) is found to be weakly dependent on pseudorapidity in central colli- 
sions; for peripheral collisions the values of V2{f]) gradually decrease as the pseudorapidity in- 
creases. The results presented here provide further input to the theoretical models of relativistic 
nucleus-nucleus collisions and will aid in determining the initial conditions of the system, the 
degree of equilibration, and the transport properties of hot QCD matter produced in heavy-ion 
collisions. 
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